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(N ■ ABSTRACT 

, In this work, two aspects of the chemical evolution of 4 He in the Galaxy are considered on the 

^"^ ■ basis of a sample of disk planetary nebulae (PN). First, an application of corrections owing to the 

contamination of 4 He from the evolution of the progenitor stars shows that the He/H abundance 
by number of atoms is reduced by 0.012 to 0.015 in average, leading to an essentially flat Hc/H 
radial distribution. Second, a determination of the helium to heavy element enrichment ratio 
using the same corrections leads to values in the range 2.8 < AY/AZ < 3.6 for Y p = 0.23 and 
2.0 < AY/AZ < 2.8 for Y p = 0.24, in good agreement with recent independent determinations 
and theoretical models. 



1. Introduction 



o 

The chemical evolution of 4 Hc in the Galaxy can be studied on the basis of abundance determi- 
nations in photoionized nebulae, namely HII regions and planetary nebulae (Peimbert, 1990, 1995; 
Pagel, 1989, 1995). Also, recent work on HII regions in blue compact galaxies has given many 
informations on this problem, especially regarding the pregalactic helium abundance Y p (Pagel 
1995, Peimbert 1995). Recent determinations of this quantity are usually in the range Y p = 0.23 
to 0.24, and in fact most discussions now concentrate on the third decimal place (Olive et al. 2000, 
Buries et al. 1999, Izotov et al. 1999). On the other hand, some aspects of the chemical evolution 
of 4 He have much greater uncertainties, particularly the question of the He/H radial gradient, 
d(He/H)/di?, and the helium to metals enrichment ratio, AY/AZ. 

Planetary nebulae (PN) are also interesting in the study of these problems, since their He/H 
abundances are relatively well measured, with estimated uncertainties better than 20%, and large 
samples are available, spanning a large metallicity range. However, as the offspring of intermediate 
mass stars, these nebulae can have chemical compositions that are different from the original inter- 
stellar abundances. Therefore, in order to study the evolution of elements such as 4 He on the basis 
of these objects, it is necessary to take into account the He contamination by the PN progenitor 
stars. Previous work on these questions has introduced this contamination in an approximate way 
(Chiappini and Maciel 1994 or not at all (Maciel 1988, Maciel and Chiappini 1994). In the past 
few years, however, detailed determinations have been published of the He yields and the amount 
of this element produced during the dredge-up episodes that occur during the late evolution of 
intermediate mass stars, both as a function of the stellar mass and metallicity. 

In the present work, we make use of these calculations in order to estimate the 4 He contam- 
ination for a sample of disk planetary nebulae. Taking into account some relationships involving 
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the nebular abundances, the central star mass and the stellar mass on the main sequence, such 
contamination can be individually determined. As a consequence, the analysis of the He/H radial 
gradient and of the enrichment ratio can be performed on more solid grounds. In section 2 the 
question of the abundance gradients is considered, and in section 3 the AY/AZ ratio is discussed. 
In section 4 a simple model to account for the contamination of the progenitor stars is proposed, 
and the results are presented and discussed in section 5. 

2. Radial abundance gradients 

Radial abundance gradients of element ratios such as O/H and S/H arc well known, and amount 
to d\og(0/R)/dR ~ -0.07 dex/kpc (Macicl, 1997, 2000). These gradients can be derived from 
a variety of objects, comprising HII regions, planetary nebulae, B stars, etc. The study of these 
gradients includes a determination of their average magnitude, the existence of spatial variations 
along the disk and the time variation of the gradients, as the objects mentioned have different 
ages, spanning a few Gyr. Several mechanisms have been proposed to account for the gradients, 
and the reader is referred to Matteucci and Chiappini (1999) for a recent review on this topic. 

The existence of a He/H gradient is much more uncertain. Shaver et al. (1983) obtained a 
zero gradient for the ratio He + /H + from radio and optical measurements in galactic HII regions, 
and suggested a small gradient, d\og(Hc/H)/dR = —0.001 dex/kpc. Their derived dispersion 
was ad — 0.08, higher than the estimated observational untertainty, a b s ~ 0.01 to 0.02. Earlier 
work suggested a negative gradient, d(Hc/H)/dR ~ —0.005 to —0.008 kpc -1 , corresponding ap- 
proximately to dlog(He/H)/di? ~ —0.02 to —0.03 dex/kpc (D'Odorico et al. 1976, Pcimbcrt and 
Serrano 1980). 

Faundez-Abans and Maciel (1986) obtained an average gradient of g?(Hc/H) /dR ~ —0.005 
kpc -1 or <ilog(Hc/H)/(ii? ~ —0.02 dex/kpc for a sample of galactic PN. Pasquali and Perinotto 
(1993) obtained negligible gradients, d(He/H)/dR ~ —0.001 kpc -1 for a sample of disk PN, reach- 
ing a maximum value d(Hc/H)/dR ~ —0.004 kpc -1 only if objects close to the galactic disk (z ~ 0) 
were included, but in this case the galactocentric distance variation was small and the result was 
uncertain. Considering PN of types I and II (Pcimbcrt 1978) with distances from the galactic plane 
z < 300 pc, they obtained a gradient d(He/H)/di? ~ —0.003 kpc -1 , corresponding approximately 
to <i(log Hc/H) /dR ~ —0.009 dex/kpc. Amnucl (1993) obtained gradients generally in the range 
rflog(Hc/H)/c?i? ~ —0.002 to —0.026 dex/kpc for the different classes of his adopted classification 
system. A discussion of these results has been given by Esteban and Peimbert (1995). Maciel and 
Chiappini (1994) also obtained very small gradients, d(He/U)/dR ~ -0.0004 ± 0.0006 kpc" 1 for a 
sample of type II PN with known distances. Recently, Esteban et al. (1999) derived a very small 
gradient, dlog(He/H)/d.R ~ -0.004 ± 0.005 dex/kpc, corresponding to d(Re/U)/dR ~ -0.001 
kpc^ 1 based on He recombination lines of three HII regions (M8, M17 and Orion). This result is 
apparently independent of temperature fluctuations in the nebulae. 

All these determinations involving PN until now have not taken into account the He contami- 
nation by the progenitor stars during their evolution. It is well known that main sequence stars in 
the mass interval 1 < M/M Q < 8 have an important contribution to the He production, especially 
during the first and second dredge-up episodes (cf. Peimbert 1990). Therefore, a realistic determi- 
nation of the He/H gradient from planetary nebulae should take into account the He production 
of their progenitor stars. 

3. The helium to metals enrichment ratio 

The helium to metals enrichment ratio is an important parameter in the study of the chemical 
evolution of the Galaxy, and is usually determined from photoionized nebulae such as HII regions 
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and HII galaxies (Lequcux et al. 1979, Pcimbert 1990, 1995, Pagel et al. 1992, Izotov et al. 
1997, Thuan and Izotov 1998, Esteban et al. 1999). Recent work has also taken into account the 
fine structure in the main sequence of nearby stars (Pagel and Portinari 1998, H0g et al. 1998) 
and Hipparcos parallaxes for open clusters have been used to stress local variations of this ratio 
(Efrcmov et al. 1997). Results obtained from these methods are in the range 2 < AY/AZ < 6, 
and lower values are generally in better agreement with theoretical calculations. 

Planetary nebulae can also be used to determine the AY/AZ ratio, and earlier determinations 
of this quantity include D'Odorico et al. (1976) and Peimbert and Serrano (1980), with results in 
the range 2.2 < AY/AZ < 3.6. 

Maciel (1988) determined the pregalactic abundance and the helium to metals enrichment 
ratio using a sample of type II PN, assuming that the He contamination from the central star 
was negligible, with the result AY/AZ ~ 3.5 ± 0.3. Chiappini and Maciel (1994) made a first 
attempt to take into account the contamination from the PN progenitor stars. For these objects, 
the observed He abundances by mass Y can be written as 

Y = Yp + Xz z + AYs (1) 

(Peimbert and Torres-Peimbert, 1974, 1976), where Y p is the pregalactic abundance, Z is the 
heavy element abundance by mass and AY S is the contamination from the progenitor star. In the 
work of Chiappini and Maciel (1994), the contamination AY S was included as a free parameter 
based on calculations by Boothroyd (private communication). The adopted values were AV S = 
0.0,0.008,0.015 and 0.022, which were applied to all nebulae in the sample. From direct Y(Z) 
correlations, the derived ratio was in the range 3.4 < AY/AZ < 5.6 [see Table 3 of Chiappini and 
Maciel (1994)]. The results by Chiappini and Maciel (1994) were based on average corrections, 
applied to all nebulae, so that any differences in the individual behaviour of the nebulae were lost. 
Since there is now some information on the masses of the PN central stars and the corresponding 
masses on the main sequence, it is interesting to revise these results in order to account for the 
contamination of the progenitor stars individually. 



4. He production in intermediate mass stars 

Since the pioneering work of Rcnzini and Voli (1981), several investigators have presented detailed 
calculations of yields of the different elements in intermediate mass stars, based on recent evolu- 
tionary models (see for example Maeder 1992, van den Hoek and Groenewegen 1997, Boothroyd 
and Sackmann 1999, Marigo 2000). Regarding He abundances, these models take into account 
the dredge-up episodes that occur in intermediate mass stars, particularly the first and second 
processes. In this work, we have adopted the recent yields by van den Hoek and Groenewegen 
(1997) for stars with masses in the range 0.9 to 7M Q on the main sequence. Their results include 
all dredge up episodes and are based on evolutionary tracks of the Geneva group up to the AGB 
branch, coupled with a synthetical AGB model with hot bottom burning. The yields depend on 
the stellar mass and metallicity, and the standard model takes as parameters the mass loss scaling 
parameter, t]agb — 4, the minimum core mass for dredge up, M™ m = 0.58 and the third dredge 
up efficiency, A = 0.75 (see van den Hoek and Groenewegen 1997 for a detailed discussion). 

As a comparison with the van den Hoek and Groenewegen (1997) data, we have used the 
computations by Boothroyd (private communication, see also Boothroyd and Sackmann 1999, 
Sackmann and Boothroyd 1999) based on earlier theoretical models by Boothroyd and Sackmann 
(1988). These computations are restricted to the first and second dredge up processes, and show 
some decrease in the He contamination for masses roughly between one and three solar masses. As 
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a consequence, they produce smaller corrections for stellar masses under 4 M Q than the previous 
set, so that they are used for comparison purposes. 

The He excess by mass AY S can be estimated from the yields of the intermediate mass stars as a 
function of the stellar mass on the main sequence Mms and total metallicity Z. We have considered 
here Z = 0.020, which is appropriate for type II PN in the galactic disk, but different mctallicities 
do not significantly change our results. In this work, we made the simplifying assumption that 
most of the He excess produced by the progenitor stars and dredged up to the surface is mixed up 
in the outer layers and ejected into the nebulae. Therefore, the determination of the corrections is 
straightforward, and our derived contamination AY, corresponds to the largest possible correction 
for a given mass. 

In order to obtain the main sequence mass of the PN progenitor stars, we have adopted a 
two step procedure. First, we obtained the central star mass M c from the observed N/O nebular 
abundance, and then we have used an initial mass-final mass relation for intermediate mass stars 
in order to derive Mms- 

The existence of a relationship between the PN core mass and the N/O abundance has been 
discussed by a number of people, often with conflicting results (see Cazetta and Maciel 2000 
for some references). Until recently, theoretical models were unable to predict some measurable 
increase in the N/O and He/H ratios for PN with relatively massive progenitors, a situation 
that has changed with the extension of the models up to 7 solar masses on the main sequence 
with overshooting (Marigo 2000) or by the inclusion of hot bottom burning (van den Hock and 
Groenewegen 1997). A calibrated relation has been recently proposed by Cazetta and Maciel 
(2000) as 



where M c is in solar masses and N/O is the abundance of nitrogen relative to oxygen by number 
of atoms. We have adopted a = 0.689, b = 0.056 and c = 0.036 for -1.2 < log(N/0) < -0.26 
and a = 0.825, b = 0.936 and c = 1.439 for log(N/0) > -0.26. This relation is valid in the 
approximate interval —1 < log (N/O) < 0.2, and has been calibrated using synthetic models for 
AGB stars by Groenewegen and de Jong (1993) and Groenewegen et al. (1995), and NLTE model 
atmospheres by Mendez et al. (1988, 1992). Since the derived core masses are relatively large, 
M c > 0.67Mq, leading to relatively large main sequence masses, we will refer to this calibration 
as high-mass calibration. 

Alternatively, we have adopted a different calibration which leads to core masses and main 
sequence masses lower than the values quoted above. In fact, recent work on the mass distribution 
of PN central stars (Stasihska et al. 1997, Stasinska and Tylenda 1990) based on self-consistent and 
distance independent methods, suggest that more than 80% of the PN central stars have masses in 
the range M c ~ 0.55 to 0.65 M Q . These results are also in agreement with core masses derived by 
Zhang (1993) and with the mass-N/O abundance relation by Marigo et al. (1996). This calibration 
can be written as 



and replaces Eq. 2 for log(N/0) < —0.26. Since this calibration produces lower core masses, 
M c > 0.55 M and lower main sequence masses than the first calibration, we will refer to it as 
low-mass calibration. 

The average initial mass-final mass relation for the high-mass calibration was taken from the 
gravity distance work of Maciel and Cazetta (1997), and can be written as 



M c = a + b log(N/0) + c [log(N/0)] 2 



(2) 



M c = 0.7242 + 0.1742 log(N/0) 



(3) 
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M c = a + ai M MS + «2 M MS + a 3 M MS + a 4 M MS (4) 

where a = 0.5426, oi = 0.02093, a 2 = -0.01122, a 3 = 0.00447 and a 4 = -0.0003119 (see figure 1 
of Maciel and Cazetta 1997). This relation was based on models by Schonberner (1983), Blocker 
and Schonberner (1990), Groenewegen and de Jong (1993) and Weidemann and Koester (1983), 
and can be used in the whole range of intermediate star masses, namely 0.8 < Mms/M q < 7. 
For details the reader is referred to Maciel and Cazetta (1997). Used in conjunction with Eq. 2, 
this equation favours relatively large main sequence masses, Mms > 3M , in agreement with the 
NLTE model atmospheres of Mendez et al. (1988, 1992). 

For the low-mass calibration, we have adopted the following initial mass-final mass relation 

M c = 0.4877 + 0.0623 M M s (5) 

which favours lower masses, Mms > lAf©. This relation is in better agreement with the corre- 
sponding relation of Groenewegen and de Jong (1993) for Mms < 3M and with the relations by 
Schonberner (1983) and Blocker and Schonberner (1990) for higher masses. It also produces main 
sequence masses closer to the values originally proposed for type II PN by Peimbert (1978). 



5. Results and discussion 

5.1 THE Hc/H RADIAL GRADIENT 

We have applied the procedure outlined in section 4 to a sample of disk planetary nebulae used 
by Maciel and Quireza (1999) to study radial abundance gradients of O/H, Ne/H, S/H and Ar/H. 
This sample includes basically the objects of Maciel and Koppcn (1994), Maciel and Chiappini 
(1994), Costa et al. (1996, 1997) and a few objects not included in the previous samples. The set 
of abundances adopted by Maciel and Quireza (1999) is not a homogeneous one, in the sense that 
they have been derived by different groups using different computational procedures, such as ief's 
or a different treatment of the collisional excitation of He lines. However, as previously discussed 
by Maciel and Quireza (1999) and Maciel and Koppcn (1994), the final sample includes only the 
best determined abundances, for which at least two recent measurements are available, so that 
the adopted abundances have uncertainties as low as possible. A few objects are low excitation 
nebulae, which may have some contribution of neutral He, but this paper concerns basically with 
high excitation objects, which are located at least at 4 kpc from the galactic centre. Furthermore, 
we have preferentially selected objects studied by our own group, to keep the inhomogeneity of 
the sample at a minimum level. For details on the observational and reduction procedures by the 
IAG/USP group the reader is referred to Costa et al. (1996, 1997) and references therein. 

The initial sample contained 130 PN. We have removed 11 objects without reliable Hc/H 
or N/O abundances and 16 PN with very high N/O values, thus partially avoiding the effects 
of temperature fluctuations (cf. Peimbert 1995, Gruenwald and Viegas 1998) and also the large 
and uncertain corrections needed for these objects. In fact, Gruenwald and Viegas (1998) have 
presented some empirical evidence that the higher stellar temperatures characteristic of Type I PN 
are associated with higher helium and nitrogen abundances, increasing the effects of temperature 
fluctuations. Distances and abundances with sources are given by Maciel and Quireza (1999). The 
final sample contains 103 PN and is shown in Tables I to III, which include (i) the common name of 
the nebula, (ii) the galactocentric distance (kpc) adopting R = 7.6 kpc as in Maciel and Quireza 
(1999), see also Maciel (1993), (iii) to (v) the observed abundances of He, N and O, given as He/H, 
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log(N/0) and e(0) = log(0/H)+12, and (vi) the references for the abundances, which are listed 
after the tables. 

The uncorrected He/H abundances are shown in Fig. 1 (top panel) as a function of the galac- 
tocentric distance R. The straight line shows a least squares linear fit. It can be seen that no 
gradient is present, and that the average abundance is Hc/H ~ 0.106 ± 0.003. 
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Figure 1 - Top panel: uncorrected Hc/H abundances as a function of the galactocentric distance and least squares linear 
fit. Bottom panel: The same using corrections for the low mass calibration. 



Tabic I - 


Input data for PN 






name 


R 


He/H 


log(N/0) 


e(0) 


Rei. 


NGC 1535 


8.71 


0.096 


-1.07 


8.61 


1 


NGC 2022 


9.69 


0.107 


-0.53 


8.41 


1 


NGC 2371 


9.00 


0.101 


-0.52 


8.63 


1 


NGC 2392 


9.44 


0.092 


-0.12 


8.50 


1 


NGC 2438 


8.61 


0.104 


-0.69 


8.84 


1 


NGC 2452 


9.14 


0.111 


-0.15 


8.61 


6 


NGC 2792 


7.94 


0.115 


-0.20 


8.75 


1 


NGC 2867 


7.54 


0.116 


-0.60 


8.71 


1 


NGC 3132 


7.64 


0.121 


-0.46 


8.95 


1,2 


NGC 3195 


6.89 


0.124 


-0.52 


8.90 


1 


NGC 3211 


7.30 


0.119 


-0.77 


8.70 


1 


NGC 3242 


7.74 


0.105 


-0.75 


8.66 


1 


NGC 3587 


7.93 


0.098 


-0.30 


8.59 


1 


NGC 3918 


6.98 


0.109 


-0.50 


8.78 


1 


NGC 5307 


6.30 


0.095 


-0.74 


8.64 


1 


NGC 5882 


6.32 


0.105 


-0.59 


8.69 


1 


NGC 6210 


6.88 


0.116 


-0.97 


8.70 


1 


NGC 6309 


5.61 


0.118 


-0.84 


8.91 


1 


NGC 6439 


3.96 


0.107 


-0.37 


8.88 


1 


NGC 6543 


7.69 


0.112 


-0.78 


8.72 


1 


NGC 6563 


5.72 


0.120 


-0.23 


8.57 


1 


NGC 6565 


6.11 


0.101 


-0.44 


8.91 


1 


NGC 6567 


6.14 


0.107 


-0.72 


8.50 


1 


NGC 6572 


6.97 


0.119 


-0.57 


8.77 


1 


NGC 6578 


5.55 


0.110 


-0.71 


8.75 


1 


NGC 6629 


6.03 


0.100 


-0.88 


8.61 


1 


NGC 6720 


7.32 


0.114 


-0.55 


8.79 


1 


NGC 6790 


6.49 


0.105 


-0.62 


8.60 


1 


NGC 6818 


6.35 


0.116 


-0.64 


8.74 


1 


NGC 6826 


7.55 


0.103 


-0.93 


8.43 


1 


NGC 6879 


6.40 


0.105 


-0.72 


8.61 


1 


NGC 6884 


7.56 


0.117 


-0.70 


8.66 


1 


NGC 6886 


6.92 


0.120 


-0.35 


8.68 


1 


NGC 6891 


6.57 


0.113 


-0.97 


8.65 


1 


NGC 6894 


7.21 


0.100 


-0.40 


8.60 


1 


NGC 6905 


6.93 


0.104 


-0.58 


8.66 


1 


NGC 7009 


7.03 


0.101 


-0.63 


8.84 


1 


NGC 7026 


7.64 


0.113 


-0.45 


8.79 


1 


NGC 7027 


7.57 


0.111 


-0.37 


8.62 


1 


NGC 7354 


7.88 


0.095 


-0.43 


8.92 


5 


NGC 7662 


7.85 


0.113 


-0.66 


8.61 


1 


IC 351 


10.36 


0.100 


-0.95 


8.48 


1 


IC 418 


8.83 


0.100 


-0.72 


8.54 


1 


IC 1297 


4.81 


0.112 


-0.47 


8.89 


1 


IC 1747 


9.41 


0.112 


-0.45 


8.75 


1 



Table II 


- Input data for PN (continued) 




name 


R 


He/H 


log(N/0) 


e(0) 


Ret. 


IC 2003 


9.83 


0.094 


-0.58 


8.62 


1 


IC 2149 


8.65 


0.099 


-0.24 


8.54 


1 


IC 2165 


9.08 


0.105 


-0.37 


8.42 


1 


T" f \ C\ A A C\ 

IC 2448 


7.35 


0.095 


-0.31 


8.59 


1 


IC 2501 


7.49 


0.093 


-0.77 


8.92 


1 


IC 2621 


7.10 


0.090 


-0.42 


8.90 


1 


IC 3568 


8.68 


0.099 


-0.87 


8.57 


1 


IC 4406 


6.44 


0.116 


-0.85 


8.75 


2,4,7 


IC 4776 


4.39 


0.090 


-1.08 


8.86 


1 


IC 5117 


7.66 


0.110 


-0.57 


8.61 


1 


IC 5217 


8.56 


0.100 


-0.68 


8.70 


1 


Cn2-1 


4.03 


0.100 


-0.34 


9.00 


1 


Fg 1 


7.12 


0.108 


-0.88 


8.45 


1 


Hl-32 


4.92 


0.110 


-1.00 


8.58 


8 


Hl-44 


4.41 


0.100 


-0.21 


8.67 


8 


Hl-56 


4.01 


0.100 


-1.10 


8.72 


8 


H2-37 


5.61 


0.110 


-0.44 


8.62 


8 


Hb 12 


8.72 


0.105 


-0.74 


8.40 


1 


He2-21 


9.96 


0.116 


-0.86 


8.45 


2 


Hc2-29 


8.22 


0.108 


-0.40 


8.62 


6 


He2-37 


7.76 


0.119 


-0.59 


8.96 


1 


He2-47 


7.32 


0.110 


-0.74 


8.90 


7,9 


He2-48 


8.15 


0.106 


-0.59 


8.53 


1 


He2-55 


7.31 


0.111 


-0.56 


8.54 


1 


He2-67 


7.14 


0.120 


-0.87 


8.91 


7 


He2-99 


5.94 


0.098 


-0.96 


8.79 


1 


He2-115 


6.17 


0.100 


-0.87 


8.62 


1 


He2-118 


4.45 


0.087 


-0.52 


8.98 


9 


He2-123 


5.53 


0.117 


-0.11 


8.67 


1 


He2-158 


5.36 


0.098 


-0.97 


8.90 


1 


Hul-1 


10.70 


0.105 


-0.60 


8.68 


1 


J 320 


11.46 


0.107 


-0.87 


8.33 


1 


J 900 


9.65 


0.098 


-0.90 


8.60 


1 


K3-68 


13.59 


0.098 


-0.55 


8.11 


8 


Ml-1 


10.65 


0.114 


-0.60 


8.30 


1 


Ml-4 


9.08 


0.105 


-0.54 


8.50 


1 


Ml-5 


9.69 


0.100 


-1.00 


8.54 


1 


Ml-7 


13.20 


0.102 


-0.17 


8.71 


8 


Ml-14 


10.27 


0.099 


-0.98 


8.40 


3 


Ml-17 


12.20 


0.113 


-0.55 


8.80 


1,8 


Ml-25 


4.04 


0.117 


-0.67 


8.99 


1 


Ml-50 


3.96 


0.102 


-0.97 


8.74 


1 


Ml-54 


4.62 


0.112 


-0.19 


8.97 


1 


Ml-57 


4.87 


0.091 


-0.51 


8.96 


1 


Ml-60 


4.31 


0.105 


-0.26 


8.84 


9 
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Table III - Input data for PN (continued) 



name 



R 



He/H log(N/0) e(0) Ref. 



Ml-74 

Ml-80 

M2-2 

M2-10 

M2-27 

M3-1 

M3-4 

M3-5 

M3-6 

M3-15 

MaC 2-1 

Pel-18 

Th2-A 



6.38 

11.32 

8.81 

3.75 

5.31 

10.25 

10.62 

10.62 

8.68 

5.72 

11.36 

6.31 

6.44 



















099 
091 
100 
091 
120 
105 
123 
110 
090 
107 
111 
091 
092 












1 



•1 

•0 
•0 



60 
.28 
48 
.55 
.43 
74 
.51 
.20 
.27 
.33 
.36 
37 
60 



8.78 
8.59 
8.43 
9.00 
8.89 
8.39 
8.72 
8.29 
8.64 
8.41 
8.44 
8.92 
8.74 



1 
1 
1 
1 

8 
6 
6 
6 
1 
1 
3 
1 



1 



References: 

1 - Maciel and Koppen (1994) 

2 - Costa et al. (1996) 

3 - Costa et al. (1997) 

4 - Corradi et al. (1997) 

5 - Hajian et al. (1997) 

6 - Kingsburgh and Barlow (1994) 

7 - Pcrinotto et al. (1994) 

8 - Pcrinotto (1991) 

9 - Koppen et al. (1991). 

The introduction of the correction procedure described in section 4 does not affect these results 
appreciably. As an example, Fig. 1 (bottom panel) shows the results for the van den Hock and 
Groenewegen (1997) correction using the low mass calibration. Note that Fig. 1 gives the Hc/H 
abundances by number of atoms, as usual, so that the contamination AY S has to be converted into 
the equivalent quantity by number A (Hc/H). The different calibrations lead to similarly negligible 
gradients, which can be written as 



In fact, the main effect of the corrections is to reduce the average He/H abundances: the van den 
Hock and Gronewegen (1997) corrections lead to He/H ~ 0.091 ± 0.003 (high mass calibration) 
and He/H ~ 0.094 ± 0.003 (low mass calibration). Use of the Boothroyd and Sackmann (1999) 
data produces larger average abundances by roughly 0.006, as the corrections are generally smaller. 
Since the He contamination is a function of the stellar mass, apparently central stars with different 
masses are scattered homogeneously in the whole range of galactocentric distances, thus destroying 
any systematic variations. It can be concluded that, in view of the uncertainties involved both in the 
abundances and distances, it is unlikely that any He/H radial gradient could be presently detected 



rf(Hc/H) 
dR 



0.0000 ± 0.0004 . 



(6) 
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from planetary nebulae. The fact that some previous determinations led to a small gradient is 
probably due to the use of small samples, as in the case of the 39 PN by Faiindez-Abans and 
Maciel (1986). Determinations of the He/H gradient based on HII regions may also be affected by 
small samples. In this case, the distances are usually better known and no corrections are needed 
for stellar contamination, but the abundance determinations may depend on somewhat uncertain 
corrections for the presence of neutral helium in the nebulae (cf. Pcimbert 1979). 

The present results can be used to estimate an upper limit to the He/H gradient on the basis 
of the total dispersion <Td observed. Since ad — 0.04 in all cases considered, we have 



d(He/H) 



alR 



< — - ~ 0.004 kpc" 1 (7) 
AR 



corresponding to rflog(He/H)/di? ~ —0.02 dex/kpc for AR ~ 10 kpc. Therefore, the existence of 
a He/H radial gradient with magnitude similar to the O/H gradient is extremely unlikely, so that 
any He/H gradient should be lower than the O/H gradient by at least a factor of 3. This conclusion 
is in good agreement with the small gradients recently derived for galactic HII regions by Esteban 
et al. (1999), as discussed in section 2. 

Recent chemical evolution models support the present conclusions. Models by Matteucci 
and co-workers (Matteucci and Frangois 1989, Matteucci and Chiappini 1999) predict very small 
gradients given by d log (He/H) /dR ~ —0.0085 dex/kpc, approximately d(He/H.)/dR ~ —0.002 
kpc -1 . The recent models by Chiappini et al. (1997) imply a maximum gradient for the inner 
Galaxy of dlog(He/H)/dR = —0.006 dex/kpc, with values smaller by a factor 3 for the outer 
region, so that these results are consistent with the present investigation. 



5.2 THE He TO METALS ENRICHMENT RATIO 

We have also applied the procedure of section 4 to the PN sample of Maciel and Quireza (1999), 
taking the O/H abundance as representative of the total metallicity. We adopted the relation 
Z ~ 25 O/H, where O/H is the oxygen abundance by number relative do hydrogen, so that oxygen 
comprises 45% of the total abundances by mass. As shown by Chiappini and Maciel (1994), this 
is essentially the same as their independently derived relation for type II PN. Since it can also 
be applied to HII regions (see for example Peimbcrt 1990), the conversion from O/H to Z for all 
photoionized nebulae becomes straightforward. 

The simple linear relationship 1 is certainly valid for metal-poor objects such as blue compact 
galaxies and low metallicity HII regions, but it probably breaks up at some maximum metallicity, 
after which a more sophisticated model would be necessary. As pointed out by Chiappini and 
Maciel (1994), He abundances of PN show some tendency to flatten out for Z > 0.010, which was 
taken in their work as an upper limit for the application of Eq. 1 . Moreover, some large metallicity 
PN have larger than average N/O ratios, apparently due to ON cycling in the progenitor stars. For 
these objects the necessary corrections to the He abundance are also larger, and therefore more 
uncertain. In this work, we decided to extend a little further the limit set by Chiappini and Maciel 
(1994) and take into account all PN in our sample having metal abundances up to 10 6 O/H ~ 700, 
which corresponds approximately to the solar value, e(O) = log (O/H) + 12 = 8.83 (Grevesse and 
Sauval 1998), orZ~ 0.017 if we take Z = 25 O/H. This includes 81 objects, or about 80% of the 
PN in our sample, so that it is still a representative sample of the planetary nebulae in the galactic 
disk. 

Our main goal here is to determine the helium to heavy element enrichment ratio AY/AZ. 
The pregalactic He abundance Y p could also be obtained using the PN alone, as in some previous 
work (see for example Maciel 1988), but this quantity is clearly better determined on the basis 
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of very low mctallicity objects, as mentioned in section 3. Therefore, we have taken Y p as a 
parameter, with the values Y p = 0.23 and Y p = 0.24. For the purposes of the present investigation, 
the discussion on the third decimal place is largely irrelevant. 

We have then applied Eq. 1 to the PN sample and obtained plots of the He abundance by 
mass Y both as a function of the O/H abundance by number and the total metal abundance by 
mass Z, as shown in Figs. 2 and 3 (filled circles). Fig. 2 gives the uncorrected abundances and fits, 
and Fig. 3 shows the results using the corrections according to the van den Hoek and Groenewegen 
data. In both figures, the straight lines are least squares fits using Y p = 0.23 (dashed lines) and 
Y p = 0.24 (solid lines). 



200 



10 b O/H 
400 



600 



uncorrected abundances 

AY/ AZ = 5.60 ± 0.22 (Y = 0.23) 

AY/AZ = 4.73 ± 0.20 (Y = 0.24) 




s 

c3 



5x10 



0.01 

z 



0.015 



Figure 2 - Uncorrected He abundances by mass Y for PN (solid circles) as a function of the O/H abundances by number 
(top axis) and total heavy element abundance by mass Z (bottom axis). Also shown are the Sun (0) and HII regions 
(empty circles). The straight lines arc least squares fits for Y p — 0.23 (dashed line) and Y p — 0.24 (solid line). 



In order to have a better idea of the behaviour of the Y(Z) function at low and high metal- 
licities, Figs. 2 and 3 also include the Sun (0) and a sample of HII regions and metal poor blue 
compact galaxies (empty circles) from the compilation of Chiappini and Maciel (1994). These 
objects have not been taken into account in the determination of the linear fits, and are included 
for comparison purposes only. In fact, inclusion of these low metallicity objects would not affect 
the derived slopes by more than a few percent, since the main constraint at low mctallicities is 
defined by the pregalactic abundance Y p . 

From Figs. 2 and 3 it can be seen that the scatter of the low metallicity objects is lower 
than for the PN, which is confirmed by considering different samples of metal poor objects in 
the literature. The higher scatter in the PN data is partially due to the higher uncertainties in 
the abundances and also to the correction procedure adopted in section 4, both of which may 
be improved in the future. In fact, the determination of abundances will benefit from the use of 
more realistic photoionization models and better empirical formulae including improved ief's. The 
correction procedure can be also improved as better stellar models are available, which will allow 
a more accurate determination of the excess helium as a function of the observed abundances and 
stellar mass. The higher uncertainties in the PN abundances and the fact that the He/H gradient 
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is at most a factor of 3 lower than the O/H gradient are probably responsible for the lack of a He/H 
x R correlation, while some correlation between Y and Z (or O /H) can be observed, particularly 
from Fig. 3. In fact, the application of the correction procedure not only reduces the average He 
abundances and the derived AY/AZ ratio, but also decreases the uncertainty of the derived slopes, 
as can be seen by a comparison of Figs. 2 and 3. This effect is similar for both calibrations, being 
slightly stronger for the high mass calibration, as can be seen in the top panel of Fig. 3. 
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Figure 3 - The same as Fig. 2 using corrected abundances according to the van den Hoek and Groenewegen (1997) data 
for the high mass calibration (top panel) and low mass calibration (bottom panel). 

The detailed results of the different calibrations are shown in Table IV. For easier comparison 
with other investigations, this table also gives the slopes AF/A(0/H) and AY/AZ le relative to 
the oxygen abundance by number and by mass, respectively, which are given by 
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AY _ 1 AY _ 1 AY 

AZ 16 ~ CL45 ~AZ ~ 25 x 0.45 A(0/H) ' ^ ' 

Note that our ratio AY/AZiq is the same as the ratio AY/ AO as defined by Pcimbcrt (1995). As 
discussed by Esteban et al. (1999), the ratio AY/AZ\q is a better constraint to chemical evolution 
models, since no corrections for the remaining elements are needed. They derive ratios in the range 
AY/AZ W = 4.78 to 6.55, with an average AY/AZ W = 5.87 for Y p = 0.240 ±0.005. From Table IV 
we have AY/AZ i& ~ 5.27 for Y p = 0.24, so that the agreement is very good. 



Table IV - Results for the AY/AZ ratio 





uncorrected 


high M 


low M 


Y p = 0.23 








AY/A(0/H) 


139.9 ±5.6 


71.9 ±4.2 


89.7 ±4.7 


AY/AZ 16 


12.44 ±0.50 


6.39 ±0.37 


7.98 ±0.41 


AY/AZ 


5.60 ±0.22 


2.87 ±0.17 


3.59 ±0.19 


corr. coeff. 


0.94 


0.89 


0.91 


Y p = 0.24 








AY/A(0/H) 


118.3 ±5.1 


50.2 ±3.9 


68.1 ±4.3 


AY/AZ 16 


10.52 ±0.45 


4.47 ±0.35 


6.06 ±0.38 


AY/AZ 


4.73 ±0.20 


2.01 ±0.16 


2.73 ±0.17 


corr. coeff. 


0.93 


0.82 


0.87 



The AY/AZ ratio decreases from the range 4.7 < AY/AZ < 5.6, which is similar to the 
result of Chiappini and Maciel (1994), to the values 

AY 

2.8 < — < 3.6 {Y p = 0.23) (9) 
AY" 

2-0 < — < 2.8 {Y p = 0.24) . (10) 

These results are closer to independently derived ratios, as we have seen in section 3. If we consider 
the Boothroyd and Sackmann data, these values are increased by 35% approximately. 

Predicted values for the enrichment ratio are usually lower than AY/AZ ~ 3, a condition 
that can be fulfilled by the present results both for Y p = 0.23 and Y p = 0.24. Recent models by 
Allen et al. (1998) and also models based on two main infall episodes for the formation of the thin 
disk and halo/thick disk (Chiappini et al. 1997) predict ratios as low as AY/AZ = 1.6, which is 
closer to our lower limit for Y p = 0.24. 

It should also be noted that, while the derived AY/AZ ratios are similar to recent deter- 
minations in the literature, the average uncertainties are smaller, which is a consequence of the 
fact that a large number of objects has been included, with a larger metallicity spread, lead- 
ing to a more reliable determination. As an example, the recent value by Pagel and Portinari 
(1998) is AY/AZ = 3 ± 2, and results from blue compact dwarf galaxies by Thuan and Izotov 
(1998) are AY/AZ = 2.3 ± 1.0. An average including both calibrations of Table IV would give 
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AY/AZ = 3.2 ± 0.5 for Y p = 0.23 and AY/AZ = 2.4 ± 0.5 for Y p = 0.24. These results are 
similar to the recent determinations of Esteban et al. (1999) for the HII regions M8, M17 and 
Orion assuming temperature fluctuations characterized by a t 2 parameter greater than zero. It 
may be remarked that both the uncorrected results of Table IV and the ratios derived by Esteban 
et al. (1999) for t 2 = 0.0 are larger by roughly 50%. According to Peimbcrt (1995), temperature 
fluctuations and the possibility of grain condensation imply some increase in the oxygen abundance 
for a given value of the He abundance, decreasing the slope AY/AZ. The same effect is attained 
here by correcting the He abundance by a certain amount AY S . Tables I to III include several 
objects for which relatively large temperature fluctuation parameters t 2 arc estimated, such as 
NGC 2392, NGC 3211 and NGC 3242 (Liu and Danziger 1993). These objects are a small fraction 
of the sample, and display the same dispersion in Figs. 1 to 3 as the remaining objects, confirming 
that the adopted correction procedure is also adequate for these nebulae. 

It is interesting to note from Fig. 3 that there seems to be some continuity from the very low 
metallicity blue compact galaxies (10 6 O/H < 120) to the HII regions (10 6 O/H < 300) and disk 
PN (10 6 O/H > 120). This might seem surprising, as these objects have apparently had different 
chemical evolution histories. However, it could be argued that the chemical evolution of a system is 
basically defined by its initial mass function (IMF) and star formation history (SFH). The IMF is 
now believed to be universal, or at least to have only small local variations (see for example Padoan 
et al. 1997, Maciel and Rocha-Pinto 1998). The SFH is clearly "bursty" in blue compact galaxies, 
but it is generally assumed as constant in the Galaxy. However, recent work by Rocha-Pinto et 
al. (2000) has shown that our Galaxy presents clear evidences of past enhanced star formation 
periods, or bursts, followed by periods of depressed star formation, or lulls. As a consequence, 
only the average star formation rate can be considered as constant, so that the similarity in the 
chemical evolution of the different systems shown in Fig. 3 is probably not surprising at all. 

Acknowledgements 

This work was partially supported by CNPq and FAPESP. 

REFERENCES 

Allen, C, Carigi, L. and Peimbcrt, M.: 1998, Astrophys. J. 494, 247. 
Amnuel, P.R.: 1993, Monthly Notices Roy. Astron. Soc. 261, 263. 
Blocker, T. and Schonberner, D.: 1990, Astron. Astrophys. 240, Lll. 
Boothroyd, A.I. and Sackmann, I.-J.: 1988, Astrophys. J. 328, 653. 
Boothroyd, A.I. and Sackmann, I.-J.: 1999, Astrophys. J. 510, 232. 

Buries, S., Nollett, K.M., Truran, J.N. and Turner, M.S.: 1999, Phys. Rev. Lett. 82, 4176. 
Cazetta, J.O. and Maciel, W.J.: 2000, Rev. Mexicana Astron. Astrof. 36, 3 
Chiappini, C. and Maciel, W.J.: 1994, Astron. Astrophys. 288, 921. 
Chiappini, C, Matteucci, F. and Gratton, R.: 1997, Astrophys. J. 477, 765. 

Corradi, R.L.M., Perinotto, M., Schwarz, H.E. and Claeskens, J.F.: 1997, Astron. Astrophys. 322, 
975. 

Costa, R.D.D., Chiappini, C, Maciel, W. J. and Freitas-Pacheco, J. A.: 1996, Astron. Astrophys. 
Suppl. 116, 249. 

Costa, R.D.D., Chiappini, C, Maciel, W. J. and Freitas-Pacheco, J. A.: 1997, Advances in stellar 

evolution, ed. R.T. Rood, A. Rcnzini, Cambridge, 159. 
D'Odorico, S., Peimbert, M. and Sabbadin, F.: 1976, Astron. Astrophys. 47, 341. 
Efremov, Y.N., Schilbach, E. and Zinnecker, H.: 1997, Astron. Nachr. 6, 335. 
Esteban, C. and Peimbert, M.: 1995, Rev. Mexicana Astron. Astrof. SC 3, 133. 



15 



Esteban, C, Peimbert, M., Torres-Peimbert, S. and Garcia- Rojas, J.: 1999, Rev. Mexicana Astron. 
Astro}. 35, 65. 

Faiindez-Abans and M., Maciel, W.J.: 1986, Astron. Astrophys. 158, 228. 
Grevesse, N. and Sauval, A.J.: 1998, Space Sci. Rev. 85, 161. 
Groenewegen, M.A.T. and dc Jong, T.: 1993, Astron. Astrophys. 267, 410. 
Groenewegen, M.A.T. , van den Hoek, L.B. and de Jong, T.: 1995, Astron. Astrophys. 293, 381. 
Grucnwald, R.B. and Viegas, S.M.M.: 1998, Astrophys. J. 501, 221. 
Hajian, A.R., Balick, B., Terzian, Y. and Pcrinotto, M.: 1997, Astrophys. J. 487, 313. 
H0g, E., Pagel, B.E.J. , Portinari, L., Thejll, P.A., Macdonald, J. and Girardi, L.: 1998, Space Sci. 
Rev. 84, 115. 

Izotov, Y.I., Thuan, T.X. and Lipovetsky, V.A.: 1997, Astrophys. J. Suppl. 108, 1. 
Izotov, Y.I., Chaffee, F.H., Foltz, C.B., Green, R.F., Guseva, N.G. and Thuan, T.X.: 1999, Astro- 
phys. J. 527, 757. 

Kingsburgh, R.L. and Barlow, M.J.: 1994, Monthly Notices Roy. Astron. Soc. 271, 257. 
Koppen, J., Acker, A. and Stenholm, B.: 1991, Astron. Astrophys. 248, 197. 
Lcqucux, J., Peimbert, M., Rayo, J.F., Serrano, A. and Torres-Peimbert, S.: 1979, Astron. Astro- 
phys. 80, 155. 

Liu, X.W. and Danziger, J.: 1993, Monthly Notices Roy. Astron. Soc. 263, 256. 
Maciel, W.J.: 1988, Astron. Astrophys. 200, 178. 
Maciel, W.J.: 1993, Astrophys. Space. Sci. 206, 285. 

Maciel, W.J.: 1997, IAU Symp. 180, cd. H.J. Habing, H.J.G.L.M. Lamcrs, Kluwer, 397. 

Maciel, W.J.: 2000, Chemical evolution of the Milky Way: stars versus clusters, ed. F. Giovannelli, 

F. Matteucci, Kluwer (in press) 
Maciel, W.J. and Cazetta, J.O.: 1997, Astrophys. Space Sci. 249, 341. 
Maciel, W.J. and Chiappini, C: 1994, Astrophys. Space Sci. 219, 231. 
Maciel, W.J. and Koppen, J.: 1994, Astron. Astrophys. 282, 436. 
Maciel, W.J. and Quireza, C: 1999, Astron. Astrophys. 345, 629. 

Maciel, W.J. and Rocha-Pinto, H.J.: 1998, Monthly Notices Roy. Astron. Soc. 299, 889. 
Maeder, A.: 1992, Astron. Astrophys. 264, 105. 

Marigo, P., Bressan, A. and Chiosi, C: 1996, Astron. Astrophys. 313, 545. 

Marigo, P.: 2000, Chemical evolution of the Milky Way: stars versus clusters, ed. F. Giovannelli, 

F. Matteucci, Kluwer (in press) 
Matteucci, F. and Chiappini, C: 1999, Chemical evolution from zero to high redshift, ed. J.R. 

Walsh, M.R. Rosa, Springer, 83 
Matteucci, F. and Francois, P.: 1989, Monthly Notices Roy. Astron. Soc. 239, 885. 
Mcndcz, R.H., Kudritzki, R.P., Herrero, A., Husfeld, D. and Groth, H.G.: 1988, Astron. Astrophys. 

190, 113. 

Mendez, R.H., Kudritzki, R.P. and Herrero, A.: 1992, Astron. Astrophys. 260, 329. 
Olive, K.A., Steigman, G. and Walker, T.P.: 2000, Phys. Rep. (in press) 

Padoan, P., Nordlund, A. and Jones, B.J.T.: 1997, Monthly Notices Roy. Astron. Soc. 288, 145. 
Pagel, B.E.J.: 1989, Evolutionary phenomena in galaxies, ed. J.E. Beckman, B.E.J. Pagel, Cam- 
bridge, 368. 

Pagel, B.E.J. : 1995, The light element abundances, ed. P. Crane, Springer, 155. 

Pagel, B.E.J, and Portinari, L.: 1998, Monthly Notices Roy. Astron. Soc. 298, 747. 

Pagel, B.E.J. , Simonson, E.A., Tcrlcvich, R.J. and Edmunds, M.C.: 1992, Monthly Notices Roy. 

Astron. Soc. 255, 325. 
Pasquali, A. and Perinotto, M.: 1993, Astron. Astrophys. 280, 581. 
Peimbert, M.: 1978, IAU Symp. 76, ed. Y. Terzian, Reidel, 215. 



1G 



Peimbert, M.: 1979, IAU Symp. 84, cd. W.B. Burton, Rcidcl, 307. 
Pcimbcrt, M.: 1990, Rep. Prog. Phys. 53, 1559. 

Peimbert, M.: 1995, The light element abundances, ed. P. Crane, Springer, 165. 
Peimbert, M. and Serrano, A.: 1980, Rev. Mexicana Astron. Astrof. 5, 9. 
Peimbert, M. and Torres-Peimbert, S.: 1974, Astrophys. J. 193, 327. 
Peimbert, M. and Torres-Peimbert, S.: 1976, Astrophys. J. 203, 581. 
Pcrinotto, M.: 1991, Astrophys. J. Suppl. 76, 687. 

Perinotto, M., Purgathofcr, A., Pasquali, A. and Patriarchi, P.: 1994, Astron. Astrophys. Suppl. 
107, 495. 

Renzini, A. and Voli, M.: 1981, Astron. Astrophys. 94, 175. 

Rocha-Pinto, H.J., Scalo, J., Maciel, W.J. and Flynn, C: 2000, Astrophys. J. Lett. 531, L115. 
Sackmann, I. -J. and Boothroyd, A.I.: 1999, Astrophys. J. 510, 217. 
Schonberncr, D.: 1983, Astrophys. J. 272, 708. 

Shaver, P.A., McGee, R.X., Newton, L.M., Danks, A.C. and Pottasch, S.R.: 1983, Monthly Notices 

Roy. Astron. Soc. 204, 53. 
Stasihska, G. and Tylenda, R.: 1990, Astron. Astrophys. 240, 467. 
Stasinska, G., Gorny, S.K. and Tylenda, R.: 1997, Astron. Astrophys. 327, 736. 
Thuan, T.X. and Izotov, Y.I.: 1998, Space Sci. Rev. 84, 83. 

van den Hoek, L.B. and Groenewegen, M.A.T.: 1997, Astron. Astrophys. Suppl. 123, 305. 
Weidemann, V. and Koester, D.: 1983, Astron. Astrophys. 121, 77. 
Zhang, C.Y.: 1993, Astrophys. J. 410, 239. 



